Entanglement evolution in multipartite cavity-reservoir systems under local unitary operations 
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We analyze the entanglement evolution of two cavity piiotons being affected by the dissipation of two indi- 
vidual reservoirs. Under an arbitrary local unitary operation on the initial state, it is shown that there is only 
one parameter which changes the entanglement dynamics. For the bipartite subsystems, we show that the en- 
tanglement of the cavity photons is correlated with that of the reservoirs, although the local operation can delay 
the time at which the photon entanglement disappears and advance the time at which the reservoir entanglement 
appears. Furthermore, via a new defined four-qubit entanglement measure and two three-qubit entanglement 
measures, we study the multipartite entanglement evolution in the composite system, which allows us to ana- 
lyze quantitatively both bipartite and multipartite entanglement within a unified framework. In addition, we also 
discuss the entanglement evolution with an arbitrary initial state. 



PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.Mn 

I. INTRODUCTION 

As one of the most subtle phenomena in many-body sys- 
tems, quantum entanglement has now been an important phys- 
ical resource widely used in quantum communication and 
quantum computation filial- Therefore it is fundamental to 
characterize the entanglement nature in quantum systems, es- 
pecially at a quantitative level. Till now, although bipartite 
entanglement is well understood in many aspects, the entan- 
glement in multipartite systems is far from clear and thus de- 
serve further exploration. 

Entanglement dynamical behavior is an important property 
in practical quantum information processing. This is because 
entanglement is fragile and always decays due to unwanted 
interactions between the system and its environment. A the- 
oretical study of two-atom spontaneous emission shows that 
entanglement does not always decay in an asymptotic way and 
it can be corrupted in a finite time [3J, which is referred to as 
entanglement sudden death (ESD). Some earlier studies also 
pointed out this fact that even a very weakly dissipative envi- 
ronment can disentangle the quantum system in a finite time 
|0-[3l • The ESD phenomenon has recently received a lot of 
attentions f§S]3\ (see also a review paper flu] and references 
therein), and, experimentally, it has been detected in photon 
ifnil and atom systems ifisll . 

A deep understanding on the ESD phenomenon concerns 
the problem where the lost entanglement goes. To answer the 
question, it is proper to enlarge the system to include its en- 
vironment. Recently, Lopez et al analyzed the entanglement 
evolution in a composite system consisting of entangled cav- 
ity photons with individual reservoirs 1191] . and show that the 
entanglement sudden birth (ESB) of reservoir-reservoir sub- 
system must happen whenever the ESD of cavity-cavity sub- 
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system occurs. Moreover, in Ref. ll20ll . Bai et al presented a 
entanglement monogamy relation in multipartite systems and 
analyzed quantitatively the bipartite entanglement transfer in 
the multipartite cavity-reservoir system. 

However, in the above analysis, the multipartite entangle- 
ment in the composite cavity -reservoir system is not well char- 
acterized, although the residual entanglement [20] can indi- 
cate its existence. Moreover, the authors only consider the 
symmetric initial state like |0) = ajOO) + When the 

initial state is asymmetric, the entanglement evolution can be 
very different. For example, a ax operation acting on the sym- 
metric state can change the evolution of the entangled cavity 
photons from the ESD route to the asymptotic decay route, al- 
though the two kinds of initial states have the equal entangle- 
ment. Therefore, it is desirable to consider the entanglement 
dynamical behavior for the asymmetric case and, particularly, 
find a good entanglement measure to characterize the genuine 
multipartite entanglement evolution. 

In this paper, for the asymmetric initial state modulated by 
an arbitrary local unitary (LU) operation, we analyze its entan- 
glement evolution in the multipartite cavity-reservoir system. 
In Sec. II, we derive the effective output state under the LU 
operation, in which there is only one parameter affecting the 
entanglement dynamics. In Sec. Ill, we analyze the bipar- 
tite entanglement transfer in the composite system, and point 
out the cavity photon entanglement is still correlated with the 
reservoir entanglement although the local operation can delay 
the ESD time and advance the ESB time. In Sec. IV, the mul- 
tipartite entanglement evolution is studied via a new defined 
four-qubit entanglement measure and two three-qubit entan- 
glement measures. In Sec. V, within a unified framework, we 
investigate the relation between bipartite entanglement trans- 
fer and multipartite entanglement transition in the composite 
system. Finally, we discuss the entanglement evolution with 
an arbitrary initial state and give a brief conclusion in Sec. VI. 
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II. THE EFFECTIVE OUTPUT STATE UNDER THE LU 
OPERATION 

Before the derivation of the effective output state under the 
LU operation, we first recall the multipartite cavity-reservoir 
system. In Ref. lfl9ll . Lopez et al considered two entangled 
cavity photons being affected by the dissipation of two indi- 
vidual A^-mode reservoirs where the interaction of a single 
cavity-reservoir system is described by the Hamiltonian 

N N 

H = hwa^a + h'^^LOkb\bk + h'^^gk{ab\ + bkO^)- (1) 

k=l k=l 

The authors analyzed the entanglement evolution with the 
symmetric initial state 

(2) 

in which the reservoirs are in the vacuum state and the quan- 
tum state of cavity photons is invariant under the permutation 
of the qubits ci and C2. They show that, along the time evolu- 
tion, the BSD of two photons can happen when the initial state 
amplitudes satisfy the condition a < (3, and this procedure is 
necessarily related to the ESB of two reservoirs. 

Now, we consider the asymmetric initial state modulated by 
an arbitrary single-qubit LU operation. Without loss of gener- 
ality, we assume that the operation acts on the first cavity, and 
then the initial state can be written as 

|a>g) = (3) 

For an arbitrary single qubit LU operation, one can decom- 
pose it as [ 21] 

C/(C,?7,7,'5) = e<RMRy{l)M5). (4) 

where the e''» is a global phase shift and Rk{0) = 
exp(— j0(Tfc/2) is the rotation along the k{— y, z) axis with 
the cTfe being the Pauli matrix. In this case, the output state 
under the time evolution is 

~ Ue,rAH\t)(^Ua,rAH",t)[Ryij)cA'^0)l (5) 

Where H' = Rl{i^),,HRMc„ H" = R\{5),,HR,{5),,, 
and the ~ means the states on two sides are equivalent up to 
some LU operations (for a detail derivation, see the appendix). 
After considering the effect of the evolution UciniH' ,t) on 
the entanglement dynamics, we find that it is equivalent to 
that of the evolution UciniH ,t) (in the appendix, we give 
the proof). The case for the evolution C/c2r2 {H" , t) is similar. 
Then Eq. (5) can be rewritten as 

\^t)^U,,rAH:t)®U,,rAH.t)[Ry{l)c^\^0)]. (6) 

which means that, under an arbitrary LU operation 
Uci (Cj '7i 77 <5), the entanglement evolution is only sensitive to 
the rotation Ry{'j)c-^. 



Therefore, the effective initial state for the entanglement 
evolution is 

|*o> = i?y(7)cJ$o) = («|00)+/?|il)),,,jOO),,,,, (7) 

in which the new basic vectors are |0) = cos(7/2)|0) + 
sin(7/2)|l) and |i) = -sin(7/2)|0) + cos(7/2)|l>. For the 
output state, we use the approximation JT^ 

C/(7?,t)e.|10) =e|10>+x|01), (8) 

where the amplitudes are £^{t) = exp(— Ki/2) and x{t) = 
[1 — exp(— Ki)]^/^ in the limit of iV — > oo for a reservoir with 
a flat spectrum. Then the effective output state has the form 

\M>t) = a(cos^|OO)+sin||0t))e,.jOO),,,, 

^Pisin^m -COS^\q^t))c^rA^t)c,r,, (9) 

where \(j>t) = £,{t)\lO)cr + x(i)|01)cr the parameter 7 
being chosen in the range [0, tt]. 



III. TWO-QUBIT ENTANGLEMENT EVOLUTION UNDER 
THE LU OPERATION 

According to the effective output state \'i't) in Eq. (9), we 
can derive the density matrices of different subsystems and an- 
alyze their entanglement dynamical behaviors. We first con- 
sider the subsystem of two cavity photons, for which its den- 
sity matrix is 

PciC2{t) ^ ipi + i>2 + i>3 + i'i, (10) 

where ■0, = \ilji){tlji\ and the four non-normalized 
pure state components are — q;cos(7/2)|00) + 

asin(7/2)e|10) - /3sin(7/2),e|01) + ;9cos(7/2)e2|ll), 
IV2) - /3sin(7/2)x|00) - /3cos(7/2)Cx|10), 
iV-a) = asin(7/2)x|00) + /?cos(7/2)ex|01), and 
\ip4) — /3cos(7/2)x^|00), respectively. For the two 
reservoirs, its density matrix is similar to that of the cavity 
photons and the following relation holds 

Pri,.,(t) ^ S^^^[pc,C2{t)], (11) 

where S'j^^ exchanges the parameters ^ and x ii-^-i ^ — ^ X 
and X 0- 

Based on the previous analysis in Ref. 1I20I1 . we choose 
the square of the concurrence to characterize the two-qubit 
entanglement evolution. The concurrence is defined as ll22ll 
C(Pij) = max(0, y/Xi — y/X2 — y/Xs — aAI) with the de- 
creasing nonnegative real numbers being the eigenvalues 
of the matrix Rij = Pij{(Jy ® ^y)p*ij[<^y ® Cy). After com- 
puting the eigenvalues of the matrices i?ciC2 Rr^r2 lElt], 
we can obtain 

elegit) = 4[max(|a/3C'|-|/3^xl'cos2(7/2),0)]2, 

= 4[max(|a/3x'|-|/3?xl'cos2(7/2),0)]2<12) 
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FIG. 1: (Color online) Different two-qubit entanglement as a func- 
tion of the time evolution nt and the rotation parameter 7. 

Combining the concurrences and the expressions of ^ (t) and 
x{t) in Eq. (8), we know that, for a given asymmetric ini- 
tial state (given the parameters a, (3 and 7), the cavity pho- 
tons entanglement decreases and the reservoir entanglement 
increases along with the time evolution. 

It is still an unsolved problem whether or not the ESD of 
cavity photons and the ESB of the reservoirs are correlated in 
the asymmetric case. With the conditions CciC2(0 = ™d 
Crir2{t) — 0, we can deduce the times of the ESD and the 
ESB, which have the forms 

K a 

where the parameter k is the dissipative constant (note that 
= exp(— Kt/2) in the entanglement evolution). According 
to the two times, we can derive that the ESD of two photons 
occurs when f3 ■ cos^(7/2) > a, as is the case for the ESB 
of two reservoirs. This means that the correlation between the 
ESD and the ESB still holds for the asymmetric initial states, 
i.e., the ESB of the reservoirs must happen when the ESD of 
cavity photons occurs. 

As an example, we choose the initial state parameters as 
a — l/VTO and /? — 3/\/l0. In Fig. 1(a), the concurrence 
^ciC2 plotted as a function of the time nt and the rotation 
parameter 7. For a fixed value of nt, the photon entanglement 
increases with the parameter 7. When the 7 is given, the C^^^^ 
decreases along the time Kt. The ESD line (the purple line) is 
also plotted in the figure, where the 7 can delay the ESD time. 



It is interesting that the entanglement evolution changes to the 
asymptotical decay route before the 7 attains to the value tt, 
and the critical value is 7 = 2arccosY^l/3 « 1.91063. In Fig. 
1(b), the entanglement evolution of C^^^, is plotted, where 
the parameter 7 can increase the reservoir entanglement and 
advance the ESB time (the purple line). The critical value 
for the route transition is also 7 2arccos-\/l/3. Moreover, 
depending on the value of the 7, the ESB can manifest before, 
simultaneously and after the ESD. 

The two-qubit entanglement of subsystem ciri has the 
form 

ClrAt) = ?'x'[l + - «2)cos(7)]2. (14) 

In Fig. 1(c), the concurrence is plotted as a function of the 
parameters nt and 7. The maximum of C^j,.^(t) appears at 
the time Kt = ln2, and the entanglement decreases with the 7. 
However, the 7 does not change the entanglement of subsys- 
tem C2r2, because the rotation Ry{'^) acts on the first cavity. 
For the subsystems Cir2 and C2ri, we can get that they have 
the equal entanglement, which can be expressed as 

Cl,^{t) = 4[max(|«/3Cxl - ^x? co^^ {l / 2) ■ (15) 

In Fig. 1(d), the concurrence is plotted. When 7 = 0, the 
entanglement evolution experiences the ESD at the time nd ~ 
ln[3(3-\/5)/2] and the ESB at the time nt = ln[3(3 + V5)/2] 
(the two intersections between the purple line and the nt axis), 
then the entanglement changes asymptotically. Along with 
the increase of the parameter 7, the time window between 
the ESD and ESB decreases, and the window become a point 
when 7 = 2arccosy273 « 1.23096. After this value, both 
the ESD and the ESB phenomena disappear. 

rv. MULTIPARTITE ENTANGLEMENT EVOLUTION 
UNDER THE LU OPERATION 

Before analyzing the entanglement evolution, we first con- 
sider how to characterize the multipartite entanglement in the 
composite system. In Ref. 1 19], the multipartite concurrence 
Cn 1 24] can not characterize completely the genuine multi- 
partite entanglement, due to its nonzero value for two Bell 
states. In Ref. [20], it is shown that the genuine multipar- 
tite entanglement can be indicated by the two-qubit residual 
entanglement 

where the sum subscripts i' G {ci,ri} and j' G {c2,r2}, 
respectively. However, the entanglement monotone property 
of Mcin is not clear, even for the case of symmetric initial 
states. 

The average multipartite entanglement may quantify the 
genuine multiqubit entanglement based on much numerical 
analysis, which is defined as Il25ll 

. ^•"'"■'-f-'.-^"'""' , (17) 





FIG. 3: (Color online) (a) Qubit-block entanglement Eq-b{ci\c2r2) 
and (b) qubit-block entanglement -Eg - i, (r 1 j C2 r2 ) as a function of the 
parameters nt and 7. 



FIG. 2: (Color online) Block-block entanglement as a function of the 
parameters at and 7. 



ciriC2 can be written as 



where the = 2(1 — tr pf) is the linear entropy and the Cij 
is the concurrence. For the four-qubit cluster-class states, an 
analytical proof of the entaiiglement monotone property for 
the E,ns was given in Refs. ll26[[27ll . 

For the effective output state \'^t) in Eq. (9), we can com- 
pute the average multipartite entanglement E^s and the resid- 
ual entanglement Afciri- After comparing the two measures, 
we can obtain that they are equivalent up to a constant factor 
2. Therefore, we can define entanglement measure 



(18) 



which quantifies the genuine multipartite entanglement be- 
tween the blocks ciTi and C27'2 and its entanglement mono- 
tone property is based on the numerical analysis on the av- 
erage multipartite entanglement E,ns- In Fig. 2, the Ebb 
is plotted as a function of the parameters Kt and 7, where 
the initial state parameters are chosen as a = and 
/3 = 3/\/T0. When 7 = 0, the EBB{nt) increases from 
to 0.36 in the region nt E {0,ln(3/2)}, then it keeps invari- 
ant until the time Kt — ln3, finally, the EBBi^t) decreases 
asymptotically. With the increase of the 7, the width of the 
plateau decreases and can be expressed as 



tu, = ln[3cos2(7/2) - 1]. 



(19) 



When 7 = 2arccos(-\/2/3), the width changes to zero and the 
evolution time is Kt — ln2. After this value, the block-block 
entanglement decreases along with the 7, and vanishes when 

7 = TT. 

The genuine tripartite entanglement in the composite sys- 
tem can be quantified by the mixed state three-tangle ll28ll 



(20) 



where t(v5^^"/) = t, - Cf - Cf^ 12911 is the pure state three. 



ijk) ^ij 

tangle and the minimum runs over all the pure state decom- 
positions of pijk- The reduced density matrix of subsystem 



(21) 



where the non-normalized pure state components are 
\^pi{t)) = acos(7/2)|000) - /3sin(7/2)^|001) + 
asin(7/2)x|010) + asin(7/2)^|100) + /3cos(7/2)^2|^oi> + 
/3cosCx|011) and \^2{t)) = /3sin(7/2)x|000> - 
/3cos(7/2)x^|010) - /3cos(7/2)^x|100), respectively. 
It is obvious that the \if2) is a separable state and its three- 
tangle is zero. Moreover, for the component \ we can 
derive t((/3i) = 0. So, the decomposition in Eq. (21) is 
the optimal and the mixed state three-tangle r^ipciric^) is 
zero. Similarly, we can obtain that all the other mixed state 
three-tangles T^{pijk) are zero. 

Although all the T^{pijk) are zero in the entanglement evo- 
lution, the three-qubit states are still entangled in the qubit- 
block form isol [3TI1 . which is not equivalent to the mixed 
state three-tangle and can not be accounted for the two-qubit 
entanglement. The qubit-block entanglement characterizes 
the genuine three-qubit entanglement under bipartite cut be- 
tween a qubit and a block of qubits, and can be defined as 
Eq-BiP;\jk) = Cf^^k - - '^ik^ i" which the C,\jk quan- 
tifies bipartite entanglement between the qubits i and jk. For 
the subsystems ciC2r2 and riC2r2, their qubit-block entangle- 
ment are 



Eq-BiPci 
Eq-B{Pri 



C2r2l ^Ci\c2r2 

) - 

C2r2/ ri\c2r2 



C1C2 
,2 



C. 



0. 



cir-2 
2 



(22) 



where C~ 



4q!^;S^X^. and the 



expressions of two-qubit concurrences Cf^ are given in Eqs. 
(12), (14) and (15). 

In Fig. 3, we plot the qubit-block entanglement as a func- 
tion of the parameters Kt and 7, where the initial state param- 
eters are chosen as a = l/\/lO and /? — 3/\/l0. For a given 
value of the 7, the qubit-block entanglement E^^^^^^^ (in Fig. 
3(a)) increases first with the time kI, and then decreases with 
the Kt after attaining to its maximal value. Along with the in- 
crease of the 7, the maximal value of E^^\^^^^ decreases. For 
the qubit-block entanglement i?ri|c2r2 (in Fig. 3(b)), the trend 
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KI (c) ((j) 

FIG. 4: (Color online) Entanglement evolutions modulated by differ- 
ent parameter values: (a) 7 = 0, (b) 7 = 71, (c) 7 = 72 and (d) 
7 = TT. The entanglements Cc^ca' C'l^ira' C'cira + C'cari and Ebb 
are represented by the green dashed line, red dotted line, brown dot- 
dashed line and blue solid line, respectively. 



of entanglement evolution is similar. In Refs. ll20l [31I] . it is 
pointed out that the qubit-block entanglement comes from the 
genuine multipartite entanglement in the enlarged pure state 
system. Here, for the multipartite cavity-reservoir system, we 
can derive the following relation 

EsBi'^t) ^ Eg^B{Pci\c2r2i't)) + Eq-B{Pri\c2r2{'t))i (23) 

which means that the qubit-block entanglement comes from 
the genuine block-block entanglement in the composite sys- 
tem. 



V. ENTANGLEMENT TRANSFER AND ENTANGLEMENT 
TRANSITION UNDER THE LU OPERATION 

Because the evolution Uc^niH ,t) ® Uc2r2{H ,t) are two 
local unitary operations under the partition Ciri|c2r2, the bi- 
partite entanglement ^^^^1 1021-2 invariant, and the following 
relation holds 

^^c,n|c2.2(**) = ^BBit) + Cl^,{t) = Aa'p^ (24) 

where i' G {ci,ri} and j' e {c2,r2}, respectively. There- 
fore, in the multipartite cavity-reservoir system, we can char- 
acterize the entanglement evolution under a unified frame- 
work, where the two qubit entanglement transfer is quantified 
by the concurrence and the multipartite entanglement tran- 
sition is quantified by the block-block entanglement Ebb- 
In Fig. 4, the entanglement evolution modulated by differ- 
ent value of 7 is plotted, where 71 — 2arccos(-\/2/3), 
72 — 2arccos(^l/3), and the initial state parameters are 
a = l/VTOand/3 = 3/VW. 

In Fig. 4(a), the parameter is chosen as 7 = 0, which 
corresponds to the symmetric initial state. In the time inter- 



val Kt € (0,ln(3/2)), a part of the initial photon-photon en- 
tanglement C^j^^ first transfers to the subsystems cir2 and 
C2ri (the brown dot-dashed line where a factor 5 is multi- 
plied), then the remaining photon-photon entanglement and 
the cavity-reservoir entanglement transition completely to the 
genuine block-block entanglement Ebb (the blue solid line). 
Along with the time evolution, the block-block entanglement 
keeps invariant and is immune to the cavity-reservoir interac- 
tion in the time interval Kt e [ln(3/2), ln3]. Finally, when 
Kt e (ln3,6], the multipartite entanglement Ebb transitions 
to the two-qubit reservoir-reservoir entanglement (the red dot- 
ted line) and the cavity-reservoir entanglement. When the pa- 
rameter 7 G (0,71], the trends of bipartite and multipartite 
entanglement evolutions are similar to those when 7 = 0, 
but the immune region of the block-block entanglement de- 
creases with the parameter and the other evolution regions ex- 
tend. In Fig. 4(b), the parameter is chosen as 7 = 71, where 
the plateau region of the Ebb changes to a point (Kt = ln2) 
and the procedures of entanglement transfer and entanglement 
transition need more time. 

When the parameter 7 e (71, 72], the initial photon entan- 
glement transfers to not only the subsystems C1T2 and C2ri 
but also the subsystem riT2, and the two-qubit entanglement 
can not transition completely to the genuine block-block en- 
tanglement. As shown in Fig. 4(c), the entanglement transfer 
and entanglement transition are plotted when 7 = 72. Along 
with the increase of the 7, the decay of the photon entangle- 
ment slows down and the transfer ratio of the two-qubit en- 
tanglement increases. At the same time, the transition ratio 
of the block-block entanglement decreases. In Fig. 4(d), the 
parameter is chosen as 7 = tt, in which the transition between 
the two-qubit entanglement and the multipartite entanglement 
disappears, and the entanglement evolution consists of only 
two-qubit entanglement transfer 

It should be pointed out that, in the unified framework of 
entanglement evolution, the entanglement monotone property 
of Eb b is based on the numerical analysis on the average mul- 
tipartite entanglement Ems I25il . The analytic proof is still an 
open problem. 



VI. DISCUSSION AND CONCLUSION 

In the more general case, an arbitrary initial state has the 
forming) = (ai|00)+a2|01)+a3|10)+a4|ll))cic2|00),,,„ 
which corresponds to two LU operations [/ci (Ci i '7i 1 7i j '^1 ) ® 
Uc2 (C2, ?/2, 72, S2) acting on the symmetric initial state. In this 
case, the analytical characterization for the entanglement evo- 
lution is not available so far However, the coiTelation be- 
tween the BSD of cavity photons and the ESB of reservoirs 
still holds. This is because we can deduce the relation 

P?,c2(^>x)-^CoxK.2(^'X)], (25) 

where the evolution Ucr{H,t)\10) = $|10) + x|01) is used. 
Based on this relation, we can obtain that when the BSD 
of cavity photons occurs at the time tESD = to, the BSB 
of reservoirs will necessarily happen at the time tESB = 
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— (1/K)ln[l — exp(— Kto)]- Moreover, the entanglement evo- 
lution is restricted by the monogamy relation 

+Cir,(t) (26) 

and the multipartite entanglement can be indicated by the two- 
qubit residual entanglement Mcin [20, 32]. 

The entanglement evolution with the asymmetric initial 
state is worth to consider for other physical systems, for ex- 
ample, the atoms systems [3], quantum dots and spin chains 
etc. Ii33i435ll . Moreover, the dissipative entanglement evo- 
lution has close relation with the type of the noise environ- 
ment. Therefore, the non-Markovian environment, the cor- 
related noises, and some operator channels ll36l - [39ll are also 
worth to study in future. 

In conclusion, we have investigated the entanglement evo- 
lution of multipartite cavity-reservoir systems with the asym- 
metric initial state. It is shown that there is only one parame- 
ter in the LU operation affecting the entanglement dynamics, 
which can delay the ESD of the photons, advance the ESB of 
the reservoirs, change the evolution route of bipartite entan- 
glement, and suppress the multipartite entanglement. How- 
ever, the correlation between the ESD and the ESB still holds. 
Furthermore, by defining the block-block entanglement, we 
analyze the multipartite entanglement evolution in the com- 
posite system, which allows us to study quantitatively both the 
entanglement transfer and the entanglement transition within 
a unified framework. Finally, the entanglement evolution with 
an arbitrary initial state is discussed. 
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Appendix 

We first prove Eq. (5). The output state under the time 
evolution is 

where, in the second equation, the identity operator / = 
Uc^Ul_^ is inserted. Substituting the Uci with the expression 
Rz{ri)Ry{'-f)Rz{S) (we neglect the global phase e*'»), we can 
obtain 

= UiRl{rj)Uc,rdH,t)RMc,<»Uc,rAH,t) 

[Ryil)c^Rz{S)cA'^o)] (28) 



where the symbol ~ means the quantum states on the two 
sides are equivalent up to the local unitary operation Ui = 
[/ci^i('5)ci^J(7)ci (note that entanglement is invariant un- 
der local unitary transformation), and we use the relation 
Rl{v)c,Uc,rAH,t)Rz{v)c, = Uc,rAH\t) with H' - 
Rl{r])ciHRz{r])ci- Due to the symmetric property of initial 
state |$o), we have the relation i?z((5)ci I'I'o) = ^^2('5)c2 1'I'o)- 
Then the output state can be expressed further as 

= Uc.rAH' ,t) ^ Rz{6)c,Rl{^)c,Uc,rAH,t)Rz{d)c, 

~ Uc,rAH\t)®U,,rAH'' ,t)[Ry{l)cA^a)]. (29) 

where we insert the identity operator Rz{5)c2R\{5)c2 = 
I in the first equation and use the relation H" — 
R\{5)c2HRz{5)c2 in the second equation. 

Next, we will prove the effects of H' and H" are equiva- 
lent to that of H in the entanglement evolution. Because the 
Hilbert space of subsystem ciri is spanned by |00) , |01) , 1 10), 
the creation and annihilation operators are 

"^^(i o) ™'l«=(^o o)' ^^^^ 

respectively. With these expressions, the rotation operator 
Rz {rf) can be rewritten as 

Rziji) = A- cAa + B, (31) 

where the coefficients A — exp(i7]/2) — exp(— 177/2) and B = 
exp(— 177/2), respectively. After substituting the expression of 
Rz (rj) into the Hamiltonian H', we can derive 

H' = Rl{rj),,HRz{v)c, 

N 

= {A ■ a^a + B)^[fujcAa + H ujkblik 

k=l 

N 

+fi'^gkiabl + bka^)]{A ■ a) a + B) 

k=l 

N N 

= hojcAa + h'y^^ojkbl.bk + h^^gk{ab\ ■ e"^ 

k=l k=l 

+bka^ ■ e-'^) 
= vX{v)HVrAv), (32) 

where ^-(77) — diag{l, exp(— ^77)}, and we used the relations 
a^aa'^a = d^a = N and a^Nbk = iMl. Similarly, for the 
Hamiltonian H", we can obtain 

H" = rUS),,HRz{S),2 = vX{S)HVrAS) (33) 

with {S) — diag{l, exp(— i^)}. Therefore, the output state 
in Eq. (29) can be written as 

l*t) = VX (77)y,t {S)Uc,r^ {H, t) ® Uc.r, {H, t) 

[Ry{l)c,VrAri)VrAm^)]- (34) 
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In the above equation, the local unitary operation 
Vj!_^{r])Vj}^{6) does not change the entanglement evolu- 
tion. Moreover, due to the reservoirs being in the vacuum 
state, we have K-i (f7)K-2 ('5)|*o) = \^o)- Therefore, the 
effective output state in the entanglement evolution has the 
form 



This means that, for the asymmetric initial state modulated 
by an arbitrary LU operation 11^ (C, 7i the entanglement 
evolution is only sensitive to the rotation Ryij). 
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